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ABSTRACT 



Aims. Molecular emission lines originating in the nuclei of luminous infra-red galaxies are used to determine the physical properties 

of the nuclear ISM in these systems. 

Methods. A large observational database of molecular emission lines is compared with model predictions that include heating by UV 

and X-ray radiation, mechanical heating, and the effects of cosmic rays. 

Results. The observed line ratios and model predictions imply a separation of the observed systems into three groups: XDRs, UV- 

dominated high-density («>10''cm"') PDRs, and lower-density (n=10'*"^cm"') PDRs that are dominated by mechanical feedback. 

Conclusions. The division of the two types of PDRs follows naturally from the evolution of the star formation cycle of these sources, 

which evolves from deeply embedded young stars, resulting in high-density (n>10^cm"^) PDRs, to a stage where the gas density has 

decreased (n=IO'*^cm"^) and mechanical feedback from supernova shocks dominates the heating budget. 

Key words. ISM: molecules - Galaxies: ISM - Galaxies: starburst ~ Galaxies: active - Infrared: galaxies - Radio lines: galaxies 



1. Introduction 

Dense (n> 10'*cm"^) molecular gas plays an important role in the 
physics of (Ultra-) Luminous Infrared Galaxies. It provides the 
materials from which stars are formed and fuels possible ac- 
tive galactic nuclei. The radiation originating in molecular gas 
in the nuclei of these galaxies provides information about the 
physical properties of the nuclear environment, such as the gas 
density and chemical composition, and the dominant type of ra- 
diation field. There are two important types of radiation: UV 
radiation (6-13.6 eV) in star-forming re gions, which generates 
Photo n Dominated Regions (PDRs; e.g. iHoUenbach & TielensI 
119991) : and X-rays (1-100 keV) emanating from accreting black 
holes, which produce in X-ray Dominated Reg ions (XDRs; e.g. 
iLepp & Dalgarnolll996t lMalonevetanil996h . Analysis of the 
state of molecular gas in a significant number of galaxies allows 
systematic effects to be identified, which can provide insight into 
the processes influencing the (gas in the) nuclei of these galax- 
ies: the star- formation rate, the evolutionary stage of the system, 
an d feedback process es. 

iBaan et all ( l2008h presented data for the CO, HCN, HNC, 
HCO"^, CN, and CS line emission of 37 infrared luminous galax- 
ies and 80 additiona l sources taken from the literature (e.g., 
Nguven-0-Rieu et all [l99l ISolomon et alj 119921: lAalto et alj 



20021: iGao & Solomonll20"04l: iGracia-Carpio et alj|2006h . They 



suggested that the sources can be divided into PDRs and XDRs, 
based on the ratios of emission lines. We expand on their analy- 
sis and explore in more detail the different physical and chemical 
processes r esponsible for the diversity in line ratios found in the 
iBaan et all (|2008) data, by comparing their data with the pre- 
dictions of PD R and XDR chemistry models first presented by 
iMeiierink &~Spaans (2005) . 



2. PDR and XDR models 

Figure [T] shows the HCN, HNC and HCO^ line ratios from 
Baan et al. (2008). To identify the physical processes that gen- 
erate the line ratios measured, we compare the data with 
the predictions of a large nii mbe r of PDR and XDR mod- 
els. IMeiierink & SpaansI l2005h and lMeiierinket alj (l2007h pre- 
sented a code to model one-dimensional PDR and XDR clouds. 
It takes into account the various heating and cooling processes 
as well as the chemistry of the gas. Because the excitation of 
the molecular species is included, the intensities of individual 
molecular emission lines can be calculated as well. Our model 
predictions for both PDRs and XDRs are shown in Fig. [1] by 
solid lines: the green and cyan lines show the XDR results and 
the blue, purple, and red ones the PDRs. These lines represent 
the line ratios emanating from different depths in the cloud, and 
cover a column density range from about Nh=10 
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10-10 cm" (depending on the density of the cloud). Three 
conclusions can be drawn from this comparison: 
1: PDRs and XDRs can be separated using the HNC/HCN line 
ratio. In XDRs, HNC is always stronger than the HCN line, pro- 
ducing a line ratio higher than unity, whereas for PDRs the in- 
verse is generally true (see Fig.[TJ),c). 

2: In the case of PDR chemistry, models with different densi- 
ties can be distinguished using the ratio of HCO^ to HCN or 
HNC, which decreases with increasing density (see Fig. [T^). 
This follows from the increased dissociative recombination rate 
of HCO^with free electrons. A change in UV flux of two orders 
of magnitude produces only modest changes in the line ratios, 
because the UV field is significantly attenuated at the high col- 
umn densities where the molecules are abundant. 
3: The most significant result from this comparison is that a sig- 
nificant part of the data points has HNC/HCN and HNC/HCO+ 
line ratios that are systematically lower than those of the 
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Fig. 1. Diagnostic diagrams using integrated 
line ratios of HCN, HNC, and HCO+ versus 
each other, a) Integrated HCO^/HCN versus 
HNC/HCO+ ratios, b) Integrated HCO+/HCN 
versus HNC/HCN ratios, and c) Integrated 
HNC/HCN versus HNC/HCO+ ratios. The 
symbols X and P and the dotted lines mark 
the regions of the line ratios of XDRs and 
PDRs. Squares indicate reliable values and tri- 
angles upper and lower limits. In addition, filled 
symbols indicate OH or H2CO megamasers. 
Colored lines indicate results for FDR and 
XDR modeling, as function of column density 
(see sect. O, where the dot denotes the high- 
est column. See Table [T] and legend for de- 
tails about the models; n denotes the density in 
cm"', and Fuv and Fx, respectively, the UV and 
X-ray flux in erg s"' cm"^. 



Table 1. The important parameters for the PDR and XDR simu- 
lations and the line styles and widths used to plot the results in 
Figs.[I]and|2l 
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data points traced by the models (log(HNC/HCN)<-0.2 and 
log(HNC/HCO-")<-0.1). A lower HNC/HCO+ ratio can be ob- 
tained by decreasing the density of the gas to lO'^'^cm"^ (red lines 
in Fig.[T]i, but the HNC/HCN ratios rema in lower than can be 
explained by the PDR and XDR models. Meiie rink et all (|2006|) 
investigated the effects of the extra heating and ionization due 
to Cosmic Rays (CRs) and showed that although CRs have an 
influence on HCO^, which is sensitive to the ionization balance 
of the gas, they have little influence on the abundances of HCN 
and HNC. All this indicates that the state of the gas is influenced 
by pro cesses that were not incorporated into the lMeiierink et all 
(I2OO7I) models. 



3. Mechanical heating 

The observed low HNC/HCN ratios could be produced by 
an enhanced conversion of HNC into HCN in warm gas. At 
temperatures higher than lOOK , HNC is converted efficiently 
into H CN (HNC-hH^HCN-hH; ISchilke et alJlT99l iTalbi et"al] 
Il996l) . Dense molec ular ga s temperatures in excess of 100 K 
have been observed. lOtt et a l. (2005) presented interferometric 
maps of NH3 of the central 500 pc of the nearby starburst galaxy 



NGC253 and found clumps with kinetic temperatures of be- 
tween -150 K and - 240 K. A sini i lar stu dy of H2CO in the star- 
burst galaxy M 82 bv lMiihle et alj (l2007h found kinetic tempera- 
tures of between ~160K and ~260K and densities on the order 
of 10"* cm""*. We note that such temperatures will occur only in 
the dense gas in the centers of the galaxies in our sample. The 
large-scale low-density (n<10'*cm"-') molecular gas (as traced by 
CO) has lower temperatures. In the PDR, XDR, and CR models, 
temperatures of 100 K or higher are only reached at the edge of 
the clouds, where the abundances of HCN and HNC are too low 
to be detectable. The regions deeper inside the clouds, where the 
molecules become abundant (at Nh > 10^''^ cm"^), are cooler 
than 100 K for radiative heating (see Fig. [3]). 

Since all radiative heating processes are insufficient to heat 
gas at larger depths into the cloud, we add mechanical heating to 
the models. Sources of this type of heating are discussed in the 
following section. The mechanical heating is implemented by 
adding a depth independent heating rate (Fmech) to the total heat- 
ing budget. A new grid of models was calculated for several den- 
sities, radiation fields, and mechanical heating rates (see Table 
[Til and the results were plotted in Fig. |2l The adopted mechan- 
ical heating rates were chosen in such a way that temperatures 
deep inside the cloud (Nh>10^'^ cm"^) were between 150K and 
250 K. The mechanical heating rates are typically about 10-50 
times lower than the radiative heating rate at the edge of the 
cloud, but start to dominate the heating balance deeper inside 
the cloud (Nh>10^'^ cm"^), where they are several orders of 
magnitude higher than both the radiative and CR heating rates. 
Figure [3] illustrates how the temperature increases when a me- 
chanical heating term is added. It also demonstrates that the in- 
crease in temperature has the desired effect on the abundances of 
the species: even though the total abundances of both molecules 
in the cloud are enhanced, the abundance of HCN increases with 
respect to HNC by a factor of about three. 

The resulting line ratios are plotted in Fig. |2] and show that, 
due to a decrease in the HNC/HCN ratio, data that could not be 
explained previously are now described well by models with a 
density of 10'*^cm"^ and a range of UV fluxes and mechanical 
heating rates (red lines). Although models with higher densities 
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Fig. 2. Diagnostic diagrams that show the ef- 
fects of the added mechanical heating. See the 
legend and Table[T]for details on the individual 
models. Plot symbols are similar to Fig.[T] 
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Fig. 3. Effects of mechanical heating on the temperature of 
the gas and the abundances of HCN, HNC, and HCO+ 

(n=10'*■5cm-^ Fuv=160 ergs"' cm^^, r,„ech=3 x IQ-'^' 



erg s 



cm ). Solid lines show the results of a model with mechanical 
heating, dashed lines for a model without. 



predict line strengths that are higher then the line strengths for 
the «=10"*'^cm"^ models and should be detectable, they produce 
line ratios that are not observed (and therefore fall outside the 
range of data points plotted in Fig. |2]i. An explanation for why 
these high-density systems are not observed is that the mechani- 
cal energiesjieededjo heat the gas are unavailable (see Sect.]?]). 
iKohno et al.l (|2008|) reported strong HCN emission in the 
cores of some Seyfert galaxies ("HCN-enhanced nuclei", or 
HENs). Because of the lack of evidence for a strong nuclear star- 
burst (e.g. lack of mid-IR PAH emission), they attributed this to 
XDR chemistry. However, our results show that a weak PDR 
with mechanical heating can exhibit the same characteristics. 



4. Sources of mechanical heating 

The addition of mechanical heating to PDR clouds provides a 
satisfactory means of explaining the observations. However, we 
need to identify an available and sufficient source of energy. 
During the evolution of a starburst, there are two feedback pro- 
cesses that potentially provide sufficient mechanical energy to 



heat dense molecular gas to the required temperatures. In the 
early stage of nuclear star- formation (the first 10^-10^ years), 
YSO outflows inject mechanical energy into individual clouds. 
Then, after about lOMyr, the first massive (M*>8M0) stars be- 
come supernovae (SNe), and inject mechanical energy into the 
star-forming region, which is distributed subsequently amongst 
the PDR clouds by turbulent dissipation. Since the energy in- 
jection required per cloud is known from the chemical modeling 
(Tmech), the required YSO outflow velocities and SN rates for the 
different models can be calculated. 

To estimate the effect of YSO outflow s on the medium, we 
follow the recipe of iHeckman et al.l (1 1 9901) to determine the en- 
ergy injection rate needed to create an YSO outflow bubble as: 



dt 



SxlO^V^^v^^^nergss- 



(1) 



where rpc is the size of the bubble in pc, vgow is the outflow ve- 
locity in kms ', and n is the initial density of the medium in 
cm~^. By assuming that the outflow energy is distributed ho- 
mogeneously over a cloud, of the same size as the bubble, the 
heating rate is the energy injection rate divided by the volume 
of the cloud: Tmech = ^^^f^/^cioud- From this, the outflow veloc- 
ity needed to produce a certain Tmech in a cloud can be calcu- 
lated (using a typical cloud size of O.lpc, which corresponds to 
a column density of Nh~10~^'^ cm"^ and a extinction of Ay w 15 
magnitudes) to be 



Vflo 



(4.1 X io^^r„ 



ch" 



-y km^ 



(2) 



This yields outflow velocities of around 10 km s"' (see Table|2]i, 
which can be achieved in YSO outflows. However, this requires 
that every dense cloud has a YSO nearby, which is unlikely since 
the YSO phase is short-lived. 

To determine the SN rate required to heat clouds on a global 
scale rather than for each individual cloud, we need to consider 
the size of the total starburst region (a sphere of diameter Dsb), 
and the number of PDRs per kpc^ (/pdr)- The transfer efficiency 
(?7trans) parameterizes the dissipation of mechanical energy as 
heat as it cascades from the large scales of the SNe shock bubbles 
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Table 2. A summary of the results of the calculations of the 
energy requirements of the models, as described in Sect. |4] 
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to the individual clouds. The total energy injection rate needed 
for the entire starburst region is then: 



dE' 



SN 



dt 



= 7^SB /f 



SB 7PDR ''cloud 



Vd 



Tmechergss ' 



(3) 



By assuming a mechanical energy output of 10^' ergs per SN, 
the SN rate can then be expressed as: 



SN rate 



dEsj^/dt 
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SNy-^ 



This yields SN rates 
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(4) 

for the 
densest 



(n=10 cm 3) clouds. A more convenient measure is to relate 
the SN rate with the star-formation rate (SFR), using the initial 
mass function (IMF) (B(M) ( Dahlen & Fransson 1999) : 



SN rate 
SFR 



= A: = 






(5) 



Using a Salpeter IMF gives a value of ^ = 0.0064. To esti- 
mate a value for /pdr, we assume a region of size (Dsb) 100 
pc and an ambient density of lO-'-lO'^cm-^, which corresponds 
to a total mass for the dense molecular gas in the SB region of 
~ IO^'^-IO^'^Mq. By assuming a clumping factor of 10 and us- 
ing the same cloud size and densities as before, we obtain a 
/pdr of about 15 PDRs pc-^. We note that only the product of 
DsB and /pdr is important for the SN rates and SFRs; a larger re- 
gion with a lower PDR density will yield the same results, since 
we determine the total molecular gas mass that can be heated by 
SNe. The energy transfer efficiency T/nans is assumed to be 10%. 
A lower efficiency is not expected, since the amount of energy 
coming from the star-formation would otherwise become larger 
than the infra-red luminosities of the sources. 

Table|2]lists the results of these calculations. It shows that the 
models that trace most of the data (n - lO^'^cm-^) require mod- 
est SFRs of only around 20 Mq yr -', whereas the models with 
higher densities require high SFRs. Such high SFRs are difficult 
to achieve and will disrupt the molecular region surrounding the 
star-forming region, which may correlate with the fact that there 
are no sources observed with line ratios corresponding to me- 
chanically heated high-density PDRs. 

Both YSO outflows and SNe can provide sufficient energy to 
heat the gas, but SNe are more likely to be the dominant process 
for the sources in our sample. First of all, SNe are able to inject 
the energy on large scales, redistributing the energy amongst the 
PDRs, while YSO outflows are only able to heat single clouds. 
This requires all stars and surrounding PDR clouds to be at the 
same (evolutionary) stage. Also, the time that an ensemble of 
stars exists in the phase where SNe are produced is far longer 



than the YSO stage, which lasts only ~ 10^-10^ years, imply- 
ing that clouds heated by SN feedback are far more likely to be 
observed. 



5. Conclusions 

Although PDR, XDR or CR chemistry models can explain the 
general behavior of the HCN, HNC andJKOMine ratios ob- 
served in nearby star-forming galaxies dBaan et al.l l2008), they 
cannot explain the low HNC/HCN line ratios observed in a large 
number of systems. The addition of heating by mechanical feed- 
back provides a solution to this problem. 

Our simulations in Figs. [1] and |2] show that the observed 
sources can be split into three main groups, based on their 
HNC/HCN line ratios. First, there are a few sources that exhibit 
XDR chemistry (log(HNC/HCN)>0). The other (PDR) sources 
can be divided into two groups: a small group (composed mostly 
of the OH-MM and ULIRG systems) of sources that can be 
modeled with standard PDR chemistry and high («>10''cm-3) 
densities and a larger group of sources with lower densities 
(n-lO^'^cm--') that are heated by mechanical feedback. 

AH evidence leads to the natural conclusion that the division 
between the PDR sources is a result of the evolution of the star- 
formation cycle of a galaxy, during which the dominant heating 
process changes and density of the gas diminishes. The high- 
density (n>10^cm-3) PDRs (with HNC/HCN ratios around unity 
and weak HCO^ lines) represent the early stage of star- formation 
where the stars are forming in the densest clouds and the stellar 
UV radiation dominates the chemistry in those clouds. During 
this period most of the gas will be converted into stars and the 
increasing amount of UV radiation of the stars will gradually 
penetrate a larger volume of the (lower-density) gas. The result- 
ing lower average density leads to an increase in HCO^ with 
respect to HCN and HNC. After about lOMyr, the first mas- 
sive stars become SNe and their shocks begin to dominate the 
heating budget, leading to lower HNC/HCN line ratios. These 
two effects alter the properties of the galaxies to those of a sec- 
ond group consisting of mechanical-feedback-dominated lower- 
density PDRs. Because of the short duty-cycle of the first stage 
compared to the second, most sources in our sample are observed 
to be in the later stage of their evolution. 

It therefore appears that mechanical feedback plays a cru- 
cial role in the time dependent physical state of the dense star- 
forming molecular gas in luminous infrared galaxies. 
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